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Abstract 

The radiative heating environment for the Huygens probe near peak heating conditions for Titan 
entry is investigated in this paper. The task of calculating the radiation-coupled flowfield, 
accounting for non-Boltzmann and non-optically thin radiation, is simplified to a rapid yet 
accurate calculation. This is achieved by using the viscous-shock layer (VSL) technique for the 
stagnation-line flowfield calculation and a modified smeared rotational band (SRB) model for 
the radiation calculation. These two methods provide a computationally efficient alternative to a 
Navier-Stokes flowfield and line-by-line radiation calculation. The results of the VSL technique 
are shown to provide an excellent comparison with the Navier-Stokes results of previous studies. 
It is shown that a conventional SRB approach is inadequate for the partially optically-thick 
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conditions present in the Huygens shock-layer around the peak heating trajectory points. A 
simple modification is proposed to the SRB model that improves its accuracy in these partially 
optic ally-thick conditions. This modified approach, labeled herein as SRBC, is compared 
throughout this study with a detailed line-by-line (LBL) calculation and is shown to compare 
within 5% in all cases. The SRBC method requires many orders-of-magnitude less 
computational time than the LBL method, which makes it ideal for coupling to the flowfield. The 
application of a collisional-radiative (CR) model for detennining the population of the CN 
electronic states, which govern the radiation for Huygens entry, is discussed and applied. The 
non-local absorption tenn in the CR model is formulated in terms of an escape factor, which is 
then curve-fit with temperature. Although the curve-fit is an approximation, it is shown to 
compare well with the exact escape factor calculation, which requires a computationally 
intensive iteration procedure. 


Nomenclature 

c = velocity of light, equal to 2. 997925xl0 10 cm/s 

C, = function of rotational temperature used in the smeared band model defined in Eq. (2) 
Ej = rotational tenn energy (cm 1 ) 

Ey = vibrational term energy (cm 1 ) 

E e = electronic tenn energy (cm' 1 ) 

e = electron charge, equal to 4. 80298xl0" 10 ( cm 3/2 -g 1/2 -s 1 ) 

F v = incoming radiative intensity integrated over all directions (erg/cm ) 
g x = degeneracy of the energy mode x 
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h = Planck’s constant, equal to 6.6256x1 0"" erg • s 

I v = frequency dependent radiative intensity ( erg • cm 2 • sr 1 ) 

j = frequency integrated emission coefficient ( erg • cm"’ -s'- sr -1 ) 

j v = frequency dependent emission coefficient ( erg • cm 3 • sr~‘ ) 

J = rotational quantum number 

k = Boltzmann constant, equal to 1.38054xl0‘ 16 erg/K 

kf r = forward rate for collisional excitation defined in Eq. (6) ( cm 3 • (g - mol) -1 • s -1 ) 

kb, r = backward rate for collisional excitation defined in Eq. (7) ( cm 3 • (g - mol) -1 • s -1 ) 

m = electron mass, equal to 9.1091x10'“ g 

N x = number density of either a specie or a molecular state specified by the subscript x 
(particles/cm ) 

N a = Avogadro’s number, equal to 6.023x10“ particles/g-mole 

Q x = partition function of the energy mode x (nondimensional) 

q r = radiative heat flux, calculated in units of erg • cm' 2 • s -1 , but usually presented in units 

of W -cm' 2 , where 1 (e/'g • cm' 2 -s -1 ) = 10 -7 (JE -cm' 2 ) 
q c = convective heat flux ( W • cm' 2 ) 

Simpj = the rate of excitation of the level i due to impact with heavy particles and electrons 
(particles- cm' 3 -s -1 ) 

Smd.i = rate of excitation of level i due to radiative emission and absorption 
(particles- cm' 3 -s -1 ) 

1/2 

T a = dissociation temperature (K), equal to (T tr T ve ) 

T tr - translational-rotational temperature (K) 

3 

American Institute of Aeronautics and Astronautics 



T ve = vibrational-electronic temperature (K) 

V = vibrational quantum number 

z = distance along the stagnation line or through a constant-property slab (cm) 

Az = thickness of a constant-property slab (cm) 

A, = escape factor for electronic state i (nondimensional) 

A, = wavelength (cm) 

v = frequency (s' 1 ) 

Vv'v" - frequency of a vibrational band head obtained from Eq (B-7) of Ref. 32 (s' 1 ) 
t, = radiative lifetime for a transition with an upper state i (s) 

Subscripts 

e = refers to the electronic energy mode 

inf = refers to a free-stream value 

J = refers to the rotational energy mode 

V = refers to the vibrational energy mode 

v = indicates frequency dependence 
Superscripts 

= indicates the radiative flux or intensity directed towards the vehicle wall 
= defines a value in the upper electronic state of a transition 
= defines a value in the lower electronic state of a transition 

I. Introduction 

1 2 

RECENT studies concerning the aerothermodynamic environment of the Huygens probe ’ , for 
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entry into Titan, have indicated that significant uncertainty exists in the prediction of shock- 
layer radiative heating. A review of these studies leads to the conclusion that this uncertainty 
may be attributed to three primary sources: 1) the accuracy of the CN violet and red molecular 
band spectral representations (these band systems contribute the majority of the radiation at 
Huygens entry conditions); 2) the influence of radiation-flowfield coupling; 3) the accuracy of 
the kinetic scheme required for modeling the CN electronic state populations. The first of these 
was made clear by Hollis et al. 1 , who showed that two widely used radiation codes, NEQAIR 
and RAD/EQUIL, disagreed by a factor of two for Huygens conditions. Wright et al. have since 
shown agreement between the line-by-line approach of NEQAIR and other similar but 
independent line-by-line codes, which suggests that the lack of agreement shown by Hollis was 
due to inadequacies in the molecular band modeling of the RAD/EQUIL code. The uncertainty 
in the spectral modeling is therefore small if the line-by-line approach is used. The drawback to 
this approach is that it is very computationally expensive, which makes it very difficult to apply 
to a coupled radiation-flowfield analysis. Therefore, although it models the spectrum accurately 
for given flowfield conditions, the inaccuracy due to ignoring coupling remains (which was the 
second source of uncertainty listed above). The coupling effect has been approximately treated in 
these past studies using a correction factor 3 developed for Jupiter entry. Unfortunately, this 
approximate method has not been validated for Titan entry with non-optically thin radiation. 
Further complications to the radiation calculation arise from the recent suggestion 4 that the 
excited electronic states of the CN molecule, which are required for the radiation calculation, are 
not populated by a Boltzmann distribution. Removing the Boltzmann assumption requires a 
kinetic scheme, similar to that used for nonequilibrium chemistry modeling, to calculate the CN 
excited state populations. Kinetics schemes for CN excitation have been proposed recently by 
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both Raynaud et al . 5 and Magin et al . 6 specifically for Titan. It has been shown that 
implementing these models reduces the radiative heating significantly 4 ' 6 . The influence of 
radiative absorption on the excitation calculation, as discussed by Bose et al. 4 , complicates the 
calculation further, especially if the spectrum is modeled using the line-by-line approach. 

The aim of the present study is to formulate the Huygens probe radiative heating prediction, 
including all effects mentioned above, into a manageable calculation. A thermo-chemical 
nonequilibrium viscous-shock-layer analysis will be used for the flowfield calculation. The line- 
by-line method of calculating the radiative spectrum will be compared with less computationally 
intensive calculation methods at conditions of interest for the Huygens probe. A new approach 
to the spectrum calculation, which is a slight modification of the smeared-rotational band 
approach, will be discussed and validated. The radiative heating to the Huygens probe near peak 
heating conditions will be calculated using the developed flowfield and radiative models for both 
the coupled and uncoupled radiation cases. Finally, a method of implementing a collisional- 
radiative model to calculate the population of the CN electronic states will be presented and its 
effect on the radiative heating at Huygens peak heating conditions will be discussed. 

II. The Viscous Shock-Layer Stagnation Line Flowfield 

The flowfield was modeled in the present study using the viscous shock layer (VSL) method 7 . 
This method allows for a computationally efficient solution that retains all of the important 
flowfield physics for the stagnation region of a blunt body. Only the stagnation line flowfield 
was considered in the present study, although the same method may be used to obtain solutions 
downstream of the stagnation line. Chemical nonequilibrium was treated following the work of 
Moss and applying the reaction rates presented by Gocken . Ionization effects have been shown 
to be negligible for Huygens entry, so the ions and electrons present in Gocken’ s model were not 
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considered. This results in a model consisting of 14 species (Ar, C, CH, CH 2 , CH 3 , CH 4 , CN, C 2 , 
H, HCN, H 2 , N, NH, N 2 ). Park’s two temperature model 10 , derived for purely nitrogen free- 
stream, was used to model vibrational nonequilibrium. This model assumes that the translational 
and rotational temperatures remain identical (T t ,) and that the vibrational and electronic 
temperatures remain identical ( T ve ). Vibrational-dissociation coupling is treated in this model by 
defining a dissociation temperature, T a = (T tr T ve ) m . The energy mode exchange and source terms 
required for the energy equations were taken from Gnoffo et al. 11 . This implementation of the 
two-temperature model into the VSL scheme followed the work of Gupta and Lee and Gupta . 
Multicomponent diffusion was modeled using the “approximate-corrected” approach presented 
by Sutton and Gnoffo 14 . The Titan atmosphere was assumed to consist of N 2 , CH 4 , and Ar in 
mole fractions of 0.9699, 0.0230, and 0.0071, respectively. An equivalent nose radius of 1.35 m 
was applied, which is slightly larger than probe’s physical nose radius of 1.25 m. This value was 
shown to give good agreement with the shock standoff and convective heating predicted by 
previous analyses 1 . A constant wall temperature of 2000 K was assumed for all cases and a wall 
emissivity of 0.9 was used for the radiation coupling calculation. This wall temperature was 
chosen to approximately match the radiative equilibrium wall values obtained in previous 
studies 1 . Therefore, the effect of radiation on the equilibrium wall temperature, which should 
have a negligible influence on the flowfield calculation, is not accounted for in this study. The 
wall was also assumed to be super-catalytic, meaning the shock-layer species at the wall were set 
equal to their free-stream values. 

The present study will concentrate on the three peak heating trajectory points identified in 
previous studies 1 . These conditions are shown in Table 1. Other trajectory points have been 
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considered throughout the course of this work, but to allow a focused and thorough discussion of 
the most severe conditions, this paper will focus exclusively on these three peak heating points. 

The flowfield solutions obtained with the VSL method developed for this study were found to 
compare well with the Navier-Stokes solutions of previous studies 1 . Figure 1 presents the 
temperature profiles for the three trajectory points. The discrete-shock used in the VSL method 
results in the translational-rotational temperature (T tr ) directly behind the shock being larger than 
what would be found using a shock-capturing method. This temperature overshoot has a 
negligible effect on both the radiative and convective heating. Figure 2 presents the CN number 
densities for the three trajectory points. 

III. Spectral Modeling of CN Molecular Band Radiation 

Various methods are available for the general modeling of molecular band spectra. The 
method chosen for determining the radiative flux from a hypersonic shock-layer is required 
primarily to model the frequency-integrated flux to the vehicle wall. Spectral details that do not 
contribute noticeably to the integrated flux are not of interest, and should be ignored if they 
require significant computational effort. The two methods that will be considered here are the 

line-by-line (LBL) approach, as presented by Arnold and Whiting 15 and used in NEQAIR 16 ' 17 , 

1 8 

and the smeared-rotational band (SRB) approach, as derived by Patch, Shackleford, and Penner 
and extended to higher order accuracy and multiple temperatures by Chambers 19,20 . The details 
of the LBL approach for a £ - £ transition, which is consistent with the CN violet (X“£ - B“£ ) 
transition, are presented in Appendix A of Ref. 32. The details of the SRB approach are 
presented in Appendix B of Ref. 32. The spectroscopic constants and absorption oscillator 
strengths (fyyj for the CN violet and red transitions were taken from Laux" for the present 
study. 
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The emission coefficients for the CN red band, and the resulting radiative flux, presented in 
this paper are Vi of the values presented in an earlier version of this paper . This is a result of the 
missing g e -lg e ■ factor in Eq. (B-2) of Ref. 32 (fortunately, this factor is equal to one for CN 
violet). The corrected values for the CN red radiation have been included throughout the present 
paper, including the radiation-coupled cases, although because of the minor importance of the 
CN red radiation for Huygens entry, this did not change any of the conclusions reached 
previously in Ref. 32. 

Assuming that it is implemented correctly and that accurate spectroscopic data is used, the 
LBL approach is considered to provide the theoretically exact radiative spectrum for a molecule. 
As a result, it provides an accurate value for the spectrally-integrated radiative flux, regardless of 
the optical depth. The drawback to this method is that it requires an enonnous number of spectral 
points (on the order of 100,000 for the spectral range of the CN violet band). This makes it very 
difficult to apply in a coupled radiation-flowfield calculation, where many radiation calculations 
are required. The SRB method, on the other hand, produces an approximate spectrum that 
captures the average shape of each vibrational band (V’ - V”). This provides a rough estimate of 
the spectrum while maintaining an accurate spectrally-integrated radiative flux in the optically 
thin limit. The accuracy of the SRB method in the optically-thin limit is a result of the method’s 
derivation, which forces this to be true for any molecular band system. Figure 3 compares the 
LBL and SRB values for the integrated emission coefficient, which is proportional to the 
integrated flux in the optically thin limit, for the CN violet band over a range of temperatures. 
The contribution from each of the major three vibrational bands is shown separately along with 
the total. In all cases, good agreement is seen between the SRB and the LBL values. This 
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confirms the statement that the SRB method is a good approximation of the LBL method for an 
optically thin gas. 

To examine how the methods compare as the gas becomes more optically thick, Table 2 
compares the integrated intensity from a 1-cm slab of CN in equilibrium at 7000 K in spectral 
range of 3. 176 to 3.444 eV for various pressures. The entire AV = 0 band and a small part of the 
AV — 1 band are captured in this spectral range. This comparison was also made by Wright et 
al. 2 , whose values are shown in the “NEQAIR96” column of the table. The NEQAIR96 values, 
which are from a line-by-line calculation, compare well with the current LBL results. A 
comparison was also made between the present LBL results and the NEQAIR96 results for the 
intensity resulting from spectral intervals of 25 Angstrom. The two values were shown to 
compare within 5% for each interval. This excellent comparison provides confidence that the 
present LBL calculation has been implemented correctly. As for the comparison between the 
LBL and SRB values, Table 2 shows that they do not compare well for the two higher pressure 
cases. This suggests that the SRB method is inadequate for modeling a gas with optically thick 
spectral regions. The optical thickness of a constant property slab is dependent on the term 
NcnAz, where Ncn is the CN number density and Az is the thickness of the slab. This tenn is a 
similarity parameter for the radiative flux emitted from a constant property slab of CN in 
equilibrium. Taking the shock-layer thickness to be about 10-cm from Figure 1 and Ncn to be 
about 1x10 cm" from Figure 2, a constant property slab representative of Huygens peak 
heating conditions is found to have a NcnAz equal to about 1x10 cm" . Figure 4 compares the 
LBL and SRB result in a small region of the spectrum (3.19-3.21 eV) for the 0.001 and 0.01 
atmosphere cases, which correspond to values of NcnAz equal to lx 10 15 and lx 10 16 cm' 2 , 
respectively. It is seen in these comparisons that the SRB model does not capture the intense 
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peaks predicted by the LBL model at each vibrational band head. This has a small influence on 
the integrated flux for the 0.001 atm case (NcnAz =1x 10 cm'") because even the largest peaks 
of the LBL spectrum are less than 75% of the blackbody limit (which is represented by nB v in 
the figure) so that the entire spectrum is essentially optically thin. For the 0.01 atm case (NcnAz 
= 1x10 cm' ), it is seen that the vibrational band heads predicted by the LBL method are 
suppressed by the blackbody limit, and much of the nearby spectrum is near this limit. These 
spectral regions are therefore optically thick. The SRB prediction for this case, on the other hand, 
fails to predict any spectral region that is more than half of the blackbody limit. The predicted 
spectrum therefore remains optically thin, which means that the absorption and corresponding 
reduction in flux predicted by the optically thick LBL method is not properly modeled by the 
SRB method. As a result, the spectrally integrated flux listed in Table 2 for the SRB model is 
larger than the LBL prediction, and the comparison between the two methods gets worse as the 
optical thickness increases. 

The previous discussion indicated that the SRB method does not provide a satisfactory model 
for the spectrum at Huygens peak heating conditions. Through the course of the present work, a 
simple correction to this model was developed that considerably improves its modeling 
capability. To rationalize the correction developed for the SRB method, it is first observed from 
the equations presented in Appendix B of Ref. 32 that decreasing the rotational temperature ( T r ) 
increases the magnitude of each vibrational band head. This implies that the spectrum has larger 
spikes, which is the behavior required for a better agreement with the LBL results. Since this 
method is an approximation to begin with, it is justifiable to consider the possibility of choosing 
a nonphysical value of T r to improve the method’s agreement with the LBL results. This is 
equivalent to introducing an empirically derived correction factor into the SRB equations 
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wherever T r appears, and does not imply that the physical value of the rotational temperature is 
believed to be any different than what is predicted by the flowfield solution (where it is assumed 
equal to the translational temperature). The initial concern in attempting this approach is that it 
may reduce the accuracy of the SRB method in the regime in which it is already known to be 
accurate, which was shown previously in Figure 3 to be the optically thin limit. Therefore, for 
this approach to be valid, it must be shown that the rotational temperature has a small influence 
on the frequency-integrated emission coefficient (j). The integrated emission coefficient from a 
single vibrational band (jyy) may be evaluated analytically for the SRB method as follows 
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and Vyv” is the frequency of the vibrational band head. The other parameters in this equation are 
defined in the Nomenclature section. The useful result of this equation is that the underlined 
tenns are the only place where the rotational temperature appears. This is significant because the 
underlined terms are negligible relative to the Vy V - term for any practical case. This is because 
Vw” is usually much larger that 1/C,- ( Vyy ~ 25,000 sec' 1 for the CN violet band while 1/C r ~ 
170 sec' 1 ). It is therefore shown that the rotational temperature has a very small influence on the 
integrated emission, which was the desired result. The next step is to find the value for T r which 
provides the best fit to the LBL calculation over a wide range of conditions. By investigating 
numerous constant property slabs over a range of temperatures and optical depths, the best value 
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was found to be T r = 1700 K. The results obtained using this approach will be referred to 
throughout this paper as the smeared rotational band corrected model, or SRBC. An example of 
the spectrum produced by this approach is shown in Figure 4. As intended, the vibrational band 
head peaks are much larger and more closely imitate the LBL spectrum. Figure 5 compares the 
flux from a constant-property slab at various optical depths and two different temperatures. The 
results consider the radiation from only the CN violet AV = 0 region of the spectrum, which is 
the most optically thick. The success of the SRBC method is clearly seen throughout the entire 
range of conditions. This result is characteristic of the other conditions investigated, which 
included temperatures ranging from 3000 to 9000 K. As mentioned previously, the constant 
property slab characteristic of the shock layer at Huygens peak heating conditions requires a 
value of NcnAz no larger than 1x10 cm' . Therefore, the SRBC model is more than adequate for 
the present application. 

IV. Uncoupled and Coupled Radiative Heating Assuming a Boltzmann Distribution 

The SRBC model was shown in the previous section to accurately predict the radiative flux 
for conditions of present interest. The use of the SRBC approach, rather than the LBL approach, 
results in a significant decrease in the computational burden associated with radiation coupling. 
But, to be completely sure that the SRBC model is providing accurate results, the SRBC and 
LBL approach will be used here for the uncoupled calculations. This allows the accuracy of the 
SRBC method to be monitored relative to the LBL calculations for an actual shock layer, and not 
just a constant property slab. For the coupled calculations, only the SRBC approach is used. 
However, once a converged coupled flowfield is obtained, the LBL approach is used to check the 
validity of the SRBC in the new flowfield. Note that no correction has been included in the 
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presented results to account for shock-layer curvature effects. In Ref. 2 a correction is applied by 
multiplying the radiative heat flux by 0.8. 

Table 3 presents a breakdown of the uncoupled radiation spectrum for the three trajectory 
points considered. The values in parentheses are the LBL results, while the values not in 
parentheses are the SRBC results. This convention will be used throughout the rest of this paper. 
Note that the CN red spectrum was only calculated using the SRBC approach. This table shows 
good agreement between the two methods for both the total flux value and the individual 
components. The distribution through the shock layer of each of these components, predicted by 
the SRBC approach for the t = 189 case, is shown in Figure 6. The corresponding radiative flux 
spectrum at the wall for the CN violet band is shown in Figure 7. It is seen that the different AV 
bands do not overlap, which means that they may be treated independently in the spectrum and 
transport calculation. This is also true for the CN red band as a whole, which is located around 1 
eV. The radiative heating values predicted by previous studies 1 are also listed in Table 3 for 
comparison. The RAD/EQUIL result is much larger than the current prediction is because it 
models the molecular band spectrum with a curve-fit over frequency, which is even more of an 
approximation than the SRB model. The blackbody limiting effect is therefore under predicted 
because the spectrum is too smooth, similarly to the SRB model. The agreement of the current 
predictions with the NEQAIR results is excellent, especially considering that different flowfield 
models were applied. The smaller contribution from the CN red band presented here, relative to 
those presented in the earlier version of this paper , improved the agreement with the NEQAIR 
results considerably. Note that the study by Olejniczak et al.~ , using NEQAIR, ignored the CN 
red contribution entirely. 
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As mentioned previously, radiation-coupled flowfields were obtained using the SRBC 
approach for the radiation model. This computationally efficient approach allowed coupled 
flowfields to be obtained without making any assumptions regarding the optical thickness, as 
was done in previous studies ’ . A discussion and comparison of the coupled flowfields will be 
presented here, followed by a discussion of the coupling effect on the radiative heating. Figure 8 
compare the temperature profiles for the t = 189 uncoupled and coupled cases. The expected 
trend for the coupled case of reduced temperature and shock standoff distance is seen. Figure 9 
compares the CN number density profiles for the coupled and uncoupled cases. The interesting 
result of increased CN number density with the addition of coupling is found. The explanation of 
this phenomenon requires the study of the contribution of each chemical reaction that involves 
CN to the total mass rate of formation of CN. Figure 10 plots the mass rate of formation of CN 
resulting from three of the main reactions involving CN. Although Gocken’s 9 simplified reaction 
set contains 6 reactions involving CN, the three not presented here do not contribute significantly 
to the difference between the coupled and uncoupled cases, although their contribution is 
included in the “Total” line presented in Figure 10. The reason for the increased CN number 
density for the coupled case is seen from this figure to be a result of the smaller negative 
formation values from reactions 1 and 3. Because these two reactions involve CN as a reactant, 
the lower temperatures throughout most of the shock layer for the coupled case results in a 
decreased forward reaction rate, and therefore lower rate of CN destruction. Reaction 2 has the 
opposite effect, but its influence is overshadowed by reactions 1 and 2. The large positive 
formation values near the shock reach almost exactly the same peak for the coupled and 
uncoupled case, which explains the similar number density profiles seen near the shock in Figure 
9. 
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Table 4 compares the convective heat flux at the wall for the coupled and uncoupled cases. A 
reduction of 1 1 to 15% is found for the three trajectory points. The influence of radiation on the 
convective heating is a combination of two competing mechanisms. The absorption of radiation 
in the boundary layer tends to increase the convective heating, while the emission of radiation 
out of the shock layer away from the body tends to decrease it. The mild decease in convective 
heating predicted here is consistent with the previous observation that the shock layer is optically 
thick only in the small spectral region of the CN violet vibrational band heads. Note that the 
values for the uncoupled case agree within 5% of the Navier-Stokes results presented by Hollis 1 . 

Table 5 presents the wall radiative heat flux values for the coupled case and compares them 
with the uncoupled values presented previously. Note that the values in parentheses are the LBL 
results that were calculated only after the coupled flowfield solution was obtained using the 
SRBC approach. An observed reduction in the flux values ranging from 18 to 25 percent 
indicates the significant influence of coupling on the radiation calculation. This reduction also 
indicates that although the CN number density was shown to increase with the addition of 
coupling, its effect on the radiation is overshadowed by the corresponding temperature decrease. 
This is because the emission coefficient depends exponentially on temperature, assuming a 
Boltzmann distribution, and only linearly on the CN number density. The column labeled 
“coupled-TW” presents the coupled radiative heating predicted using the Tauber-Wakefield 
approximate method' . This method depends on the uncoupled radiative heating and a constant, k, 
which has been assumed equal to 2 in previous studies 1 ' 2 . The present results show that this 
method over predicts the effect of coupling for the first two trajectory points. Figure 1 1 shows 
the difference in the wall directed radiative flux profiles with the addition of coupling. Although 
the magnitudes of the coupled results are lower than the uncoupled results, the shapes of the 
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profiles are similar. This indicates that the temperature decrease due to the addition of coupling 
does not introduce any fundamental changes in the radiation spectrum and transport, other than 
the expected decrease in emission. 


V. Collisional-Radiative Modeling for Titan Entry 

Recent shock-tube experiments by Bose et al . 4 have shown that near Huygens peak heating 
conditions, the measured CN(B) population* is an order-of-magnitude less than predicted by a 
Boltzmann distribution. Since the CN(B) state is the upper level of the CN violet transition, this 
implies that the radiation from the CN violet band is an order-of-magnitude less than that 
predicted assuming a Boltzmann distribution. Using excitation rates from previous studies" ’ , 
Bose et al . 4 formulated a simple collisional-radiative (CR) model that provided a closer 
prediction to their experimental results than the Boltzmann model. Raynaud et al . 5 and Magin et 
al . 6 have since assembled a more elaborate set of excitation rates and discussed their influence on 


the predicted CN(B) population. The present study will review the rates suggested by these past 
studies and consider the effect of using other values found in the literature. Also, the absorption 
tenn shown by Bose to be important, will be formulated in terms of an escape factor. It will be 
shown that this tenn can then be curve fit as a function of temperature and therefore 
implemented easily into the CR model. This avoids the costly iteration procedure used by Bose 
et al. in their treatment of the absorption tenn. The details of the CR model as implemented in 
the present study are discussed below. 

The conservation equation for the electronic state i of a molecule may be written as 


dt imp,i 


+ s 


rad,i 


particles 
cm 3 • s 


( 3 ) 


§ The notation that CN(B) refers to CN(B 2 X + ), CN(X) refers to CN(X 2 £ + ), and CN(A) refers to CN(A : n) will be used throughout this paper. 
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where the term on the left is the time derivative of the number density (particles/cm -s) of the 
molecules in the electronic state i, S un!l , is the rate of excitation of the level i due to impact with 
heavy particles and electrons, and S ra d,i is the rate of excitation of level i due to radiative 
emission and absorption. The S impj term is written as follows 


1 Nr 

S imp.i = T ~ ^ (fli.r ~ f, r N A(Y)N B{X) ~ ^b.r ^ A(2)^ B(2) ) 

^ a r = 1 


particles 
V cm 3 • s J 


( 4 ) 


where a pr and are the forward and backward stoichiometric coefficients for reaction r, 
respectively. This equation assumes that the excitation processes may be written in the following 
form 


A(l) + B(l)< >A(2) + B(2) (5) 

where A and B each represent a specie, which may be the same or different, and 1 and 2 
represent the initial and final electronic state of the specie, respectively. In Eq. (4), the N values 
with either an A or B subscript represent the number density of the particular state. Table 6 lists 
the excitation processes, and corresponding rates 5 ’ 6 , which were found to contribute noticeably to 
the CN(B) population. The electron impact excitation rates are not required for the weakly- 
ionized cases considered here, but they are listed because they may be important for other cases. 
Note that Magin et al. 6 found these reactions to have a significant effect at the present flow 
conditions because their flowfield contained a significantly higher electron number density than 
what is currently predicted (when the ionized species are considered in the flowfield). This over 
prediction of electrons was a result of their use of Nelson’s" chemical reaction rate set instead of 
Gocken’s 9 recent model. The influence of other possible values for the rates listed in Table 6 
will be discussed later. The forward rate coefficient is defined in terms of the listed parameters 
as 
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k f,r = A r T f n " exp (~C r /T f ) 


cm 


.3 A 


g - mole ■ s 


( 6 ) 


Considering reactions of the form of Eq. (4), the backwards rate coefficient may be written from 
detailed balancing as 


Kr =M r ») gA<l)gB(l) 

Sa(2)Sb(2) 


he 

~kT\ 


A(2) Ea(\) + Eb(Z) 


( 7 ) 


where g and E represent the degeneracy and term energy (1/cm) of the indicated state, 
respectively. 

The rate of excitation of level i due to radiative emission and absorption, S ra< tj, may be written 


as 


5 —^A 

° rad,i ~ yl i 

T: 


particles 
V cm 3 • s J 


( 8 ) 


where r, is the radiative lifetime in seconds, which are listed in Table 7 for CN, and A t is the 
escape factor (nondimensional). As mentioned previously, the escape factor accounts for the 
radiative absorption tenn shown by Bose et al. 4 to be important for Huygens conditions. If 
absorption is ignored, or if the gas is optically thin, the escape factor is equal to one. This escape 
factor is shown in Appendix C of Ref. 32 to be written as 


A i =x ~^~\~r~ F y dv 

Nj J 0 h v 


( 9 ) 


where K* V)t is the absorption coefficient, including induced emission, resulting from the 
molecular band system with an upper state i. The function F v is the incoming radiative intensity 
integrated over all directions, which may be written as 


F v -2n^I v {<f)sva.^d^ 


' erg 
V cm 2 J 


( 10 ) 
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2 

where /,, is the directional radiative intensity (erg/cnr-sr) and (j) is the angle from the vertical axis 
directed away from the wall. Note that if a cos <f) was added to the integrand, the conventional 
radiative flux equation would be obtained. The evaluation of this term, consistent with the 
tangent-slab approximation, and written in terms of approximate exponential integrals, is 
presented in Appendix D of Ref. 32. A significant property of the escape factor is that it depends 
on the radiation emitted from every point in the shock layer. Therefore if an exact solution is 
desired, it is necessary to iterate over the entire flowfield. Fortunately, as written in Eq. (9), it is 
not very dependent on this iteration procedure because of the large value of N, present in the 
denominator. After evaluating the escape factor exactly through an iteration procedure for the 
trajectory points considered in this study, it was found that it could be correlated as a function of 
only temperature as follows: 

A t - -4.6079 xlO -8 7^ +8.052x10^7^ -2.6814 9000K > T ve > 4477 K 

A t =0.0 T ve < 4477 K 

This correlation is compared with the computed values in Figure 12. Although some error is 
introduced by using the correlation, it allows for a considerable reduction in computational time 
and complexity, especially for radiation coupled solutions. A comparison between the radiative 
flux values computed using the exact values and the correlation is presented below. 

The solution to Eq. (3) may be obtained conveniently by applying the quasi-steady state 
assumption . This assumption assumes that derivative tenn may be set equal to zero, which 
transforms the differential equation into an algebraic equation. Both Raynaud et al. 5 and Magin 
et al. 6 investigated the accuracy of this assumption for Titan and concluded that it was 
sufficiently accurate. Combining Eqs. (3), (4), (8), and the set of excitation rates listed in Table 
6, the number densities of the A and B states of the CN molecule may be solved for, resulting in: 
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™cn(B) ~ N C n(X) 


+ kf, 4^e + kf,6^ n 2 (X,v=U) + ^f,7^ N 2 (A) 

\^b,2-^M + ^bA^e + ^ b,6^ N 2 (X ,r=\\) + ^bj^ N 2 (A) + ^CN(B)^ a ^ T CN(B) J 


( 12 ) 


and 


^CN(A) ~ N C N(X) 


kf,\N M + kf, 3^<? + k f 5 N Ni(Xy= \ i) 


kb,\N M + k b ^N e + k bs N N ^ Xv=lV) + A CN (a)N a / ?cn(A) 


(13) 


The ground state number density, CN(X), may be solved from the requirement that the sum of 
the various state number densities equals the total CN number density: 


NcN(X) - ™CN N CN(A) ^CN(B) 


(14) 


These equations may be set up as a set of three linear equations in three unknowns and solved 
easily. For the cases considered here, the A and X states are so close to their Boltzmann values 
that Eqs. (13) and (14) do not need to be considered. Therefore, only Eq. (12) is evaluated 
assuming Ncn(x) is equal to its Boltzmann value, along with all of the N 2 states. 

The influence of the various excitation mechanisms on the CN(B) number density profile for 
the t = 189 case are shown in Figure 13. The “all reactions” case with A = 1 (which means that 
absorption is not included) represents the baseline case. The influence of each reaction is then 
seen by removing each reaction individually from the calculation and observing the deviation of 
the predicted values of Ncn(B) from the baseline case. This is shown in Figure 13 for reactions 2, 
6, and 7, which are defined in Table 6. The most important reaction appears to be reaction 6, 
which is resonant-collisional deexcitation reaction. Also shown in this figure is the influence of 
including the absorption tenn, which implies using an escape factor less than one. The “A = 
exact” case implemented the exact escape factor calculation as described previously. The “A = 
fit” case applied Eq. (11) for the escape factor. It is seen that the “exact” and “fit” cases are in 
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good agreement. This is also true for the other two trajectory points considered in this study, 
which therefore validates the use of the curve-fit for these conditions. 

The wall radiative heat flux values obtained using the CR model are listed in Table 8. The 
results obtained by using the escape factor curve fit or ignoring absorption {A = 1) may be 
compared with exact escape factor case. It is seen that, again, good agreement is obtained using 
the curve fit. The reduction in q r due to the use of the CR model, rather than a Boltzmann 
distribution, is found to be as large as 70% at t = 185. The increase in q r due to absorption is 
found to be the largest at t = 193, where it provides an increase of about 16% above the A = 1 
case. The effect of radiation- flowfield coupling on the CR model was also considered and the 
resulting radiative heating values presented in the last column of Table 8. Considering the greatly 
reduced radiation predicted by the CR model, it is not surprising that the coupling effect is small. 
Note that the coupled results, which include the curve-fit model for the escape factor, are not 
even reduced back to the uncoupled cases with no absorption (A = 1). 

It was shown in this section that accounting for a non-Boltzmann distribution reduces the 
radiative heating by about half. This significant reduction is dependent on the rates used for the 
various collisional excitation processes, which are not very well known. As shown in Figure 13, 
reaction 6 is the most important while reactions 2 and 7 are of roughly equal importance. The 
rate for reaction 2 used in the present model, following past Titan studies, was taken from 
Zalogin et al. 25 This study conducted experiments in a 9.6% C0 2 , 0.3% N 2 and 90.1% Ar free- 
stream at 1 Torr and chose values for the excitation rates so that their computational and 
experimental results agreed. A concern in using the rate derived from this study is the large 
argon concentration present, which essentially implies that the collision partner in their 
experiment for reaction 2 was argon. The non-Boltzmann results obtained at relatively large 
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9Q 

initial pressures indicate that Ar may be inefficient, relative to N 2 , at exciting CN. FairbairrT 
found a similar result for the excitation of CN through collisions with argon. Past studies that 
have considered the excitation of CN through collisions with N 2 have suggested much larger rate 
constants 30,31 . Unfortunately, these studies have been at low temperature conditions (around 300 
K) so that the derived rate constants are not directly applicable to the present problem. It should 
be mentioned, however, that the rate for reactions 6-8 were obtained from similar low 
temperature conditions. Thus, it is worth considering the effect of replacing reaction 2 with one 
of these low temperature rates. If a value of 1.2x10 cm /mol/s is used for reaction 2, as was 
suggested by Tereshchenko (this is also close to value suggested by Provencher ), the resulting 
CN(B) populations are very close to a Boltzmann distribution for the cases considered here. This 
conclusion indicates that caution must be used when considering the large decrease in the 
radiative heating resulting from the present CR model. 

VI. Conclusions 

An approach to calculating the Huygens probe radiative heating was developed, which 
included the effects of radiation-flowfield coupling and non-Boltzmann electronic state 
populations. The ability to treat radiation-flowfield coupling, with reasonable computational 
requirements, was made possible by using a modified smeared rotational band (SRBC) 
calculation for the radiations spectrum. This approach was shown to provide frequency- 
integrated heat flux values within 5% of a line-by-line calculation over the range of conditions 
for Huygens entry. The reduction in radiative heating due to coupling was found to be about 20% 
for the three peak-heating trajectory points considered. This reduction was less than that 
estimated in previous studies using the Tauber-Wakefield correction. The convective heating was 
shown to be reduced by about 15%. The uncoupled radiative heating, obtained with the efficient 
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SRBC method, was shown to agree well with previous NEQAIR results. Details of implementing 
a collisional-radiative (CR) model for Titan entry were presented. The treatment of the radiation 
absorption term was formulated in terms of an escape factor. For the range of conditions 
considered, it was shown that this tenn could be approximated using a curve-fit as a function of 
temperature. The use of the CR model for calculating the population of the CN(B) state was 
shown to reduce the radiative heating by roughly a half, with the effect of coupling being small. 
However, caution is noted in accepting this conclusion because of the uncertainties in the rate 
data. 
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Tables 


Table 1. Huygens entry conditions considered in this study 


t( s) 

U inf (mls) 

Pinf (kg/m 3 ) 

T inf ( K) 

185 

5489 

1.83xl0' 4 

177.0 

189 

5126 

2.96xl0' 4 

176.6 

193 

4705 

3. 79x1 O' 4 

175.8 


Table 2. Intensity resulting from a 1-cm constant property slab 
of CN in equilibrium at 7000 K (W/cni 2 - sr) 


p (atm) 

N C n (cm 4 ) 

NEQAIR96 

Present - LBL 

Present - SRB 

0.001 

1.05xl0 15 

1.67 

1.60 

1.63 

0.01 

1.05xl0 16 

11.7 

11.6 

14.3 

0.1 

1.05xl0 17 

41.0 

39.6 

60.3 


Table 3. Breakdown of the radiative heat flux towards the wall (</,") at z = 0 for the 
uncoupled case (W/cm“), values not is parenthesis are the SRBC predictions and values 
in parenthesis are the LBL predictions. 


t (sec) 

CN violet 
AV= 0 

CN violet 
AV = 1 

CN violet 
AV = -1 

CN 

red 

Total 

NEQAIR 
(Ref. 1) 

RAD/EQUIL 
(Ref. 1) 

185 

46.7 (47.7) 

11.6(12.6) 

9.0 (9.5) 

4.4 

71.7 (74.2) 

72.3 

118.4 

189 

50.5 (50.7) 

13.5 (14.8) 

11.3 (11.8) 

7.5 

82.8 (84.8) 

81.5 

143.1 

193 

39.8 (39.3) 

11.3 (12.7) 

10.3 (10.7) 

10.5 

71.9 (73.2) 

72.1 



Table 4. Comparison of the convective heat flux for the 
uncoupled and coupled cases (W/cm 2 ) 


t (sec) 

q c 

uncoupled 

q c 

coupled 

% 

difference 

185 

43.9 

37.6 

14.4 

189 

44.5 

38.7 


193 

37.9 

33.7 

11.1 













Table 5. Comparison of the radiative heat flux towards the wall at z = 0 for the uncoupled 
and coupled cases (W/cm 2 ) 


t (sec) 

(Jr 

uncoupled 

(Jr 

coupled 

% difference 

(Jr 

coupled - TW 

185 

71.7(74.2) 

59.0 (60.8) 

17.7(18.1) 

51.8 (53.2) 

189 

82.8 (84.8) 

66.5 (68.3) 

19.7(19.5) 

61.3 (62.5) 

193 

71.9(73.2) 

54.0(55.0) 

25.0 (25.0) 

54.5 (55.3) 


Table 6. Excitation reactions considered in the present study 




A 

n 

C 

Tf 

Tb 

Ref. 

1 

Neutral Particle Collisional Deexcitation: 
CN(X) + M<r+ CN(A) + M 

1.5e+l 1 

0.5 

13300 

T a 

T tr 

Zagolin 

2 

CN(X) + M<->CN(B) + M 

1.8e+l 1 

0.5 

37000 

T a 

T tr 

Zagolin 

3 

Electron Impact Deexcitation: 

CN(X) + e'<->CN(A) + e 

6.0e+14 

0.5 

13300 

T 

1 ve 

T ve 

Zagolin 

4 

CN(X) + e'<-»CN(B) + e' 

6.3e+14 

0.5 

37000 

T ve 

T ve 

Zagolin 

5 

Resonant Collisional Deexcitation: 

CN(X) + N 2 (X, V =4) <-> CN(A) + N 2 (X,F=0) 

6.0e+13 

0 

0 


T tr 

Chemyi 

6 

CN(X) + N 2 (X, V = 1 1 ) ^ CN(B) + N 2 (X, V =0) 

6.0e+13 

0 

0 

- 

T tr 

Chemyi 

7 

Quenching Reaction: 

N 2 (A) + CN(X) ^ N 2 (X) + CN(B) 

4.2e+12 

0.5 

0 

T a 

Ttr 

Pintassilgo 


Table 7. Radiative lifetimes for the CN radiative transitions. 




T (s) 

Ref. 

CN(A) 

CN(B) 

-+CN (X) + hv 
-» CN(X) + hv 

1.54e-5 

6.55e-8 

Zagolin 

Zagolin 


Table 8. Comparison of the radiative heat flux towards the wall at z = 0 for the various 
escape factor models (W/cm") 


t (sec) 

(Jr 

uncoupled 
(A = exact) 

(Jr 

uncoupled 
(A = fit) 

(Jr 

uncoupled 
(A =1.0) 

(Jr 

coupled 
(A = fit) 

185 

21.0 

21.4 (22.0) 

18.2 

20.8 

189 

35.3 

34.6 (35.4) 

29.0 

32.8 

193 

40.3 

39.3 (39.9) 

31.1 

34.9 








Figures 



Figure 1. Stagnation line temperature profiles from the viscous shock layer solution at 

three peak heating trajectory points. 



Figure 2. Stagnation line CN number density profiles at three peak heating trajectory 

points. 
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Figure 4. Comparison between the various spectrum models for the radiative flux 
spectrum from a constant property slab of CN at 7000 K. The spectral range shown 
illustrates the (0-0) and (1-1) band heads. 



Figure 5. Comparison of the radiative flux from a constant property slab for different 

spectral models. 
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Figure 6. The contributions from the various radiative mechanisms to the wall directed 
radiative flux profde for the t = 189 uncoupled case 
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Figure 7. The CN violet spectrum at the wall as predicted by the SRBC method for the 

189 coupled case. 



Figure 8. Comparison between the coupled and uncoupled 
temperature profiles for the t = 189 case. 



Figure 9. Comparison between the coupled and uncoupled CN number density profiles. 



Figure 10. Comparison of the mass rate of formation of CN molecules by various 
reactions for the uncoupled and coupled cases at t = 189 s. 
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Figure 1 1 . Comparison between the coupled and uncoupled wall-directed radiative flux 

profdes. 
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Figure 12. Values of the escape factor calculated for the three trajectory points compared 

with an approximate curve fit. 
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Figure 13. The influence of various excitation mechanisms on the CN(B) number density 

for the t=T89 case. 


